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Abstract 

Surface reconstructions of Au film on Pd(llO) 
substrate are studied using a local Einstein ap- 
proximation to quasiharmonic theory with the 
Sutton-Chen interatomic potential. Tempera- 
ture dependent surface free energies for differ- 
ent coverages and surface structures are calcu- 
lated. Experimentally observed transformations 
from (lxl) to (1x2) and (1x3) structures can be 
explained in the framework of this model. Also 
conditions for Stranski-Krastanov growth mode 
are found to comply with experiments. The do- 
main of validity of the model neglecting mixing 
entropy is analyzed. 
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Reconstructions of (llO)-surfaces of d- 
metals have inspired numerous experimen- 
tal and theoretical studies. It has been 
confirmed by Low Energy Electron Diffrac- 
tion (LEED) [|] and Scanning Tunneling Mi- 
croscopy (STM) |2]|| that ^d-metals, (Rh, 
Pd and Ag) favour unreconstructed (lxl) 
structures whilst 5d-metals (Ir, Pt and Au) 
exhibit (1x2) reconstructions spontaneously. 
Reconstructions of pure metals are theoreti- 
cally fairly well understood either in terms of 
the relative contribution of the d-electrons to 
binding Q or in terms of the range of inter- 
atomic potentials []5]|§. However, conditions 
for surface reconstructions in heterogenous 
adsorbate-substrate systems have not been 
explained in simple terms, yet. 

A heterogenous system which has re- 
ceived both experimental and theoretical at- 
tention is Au on Ag (110). Rousset et. al. J7J 
have revealed with STM that at submono- 
layer coverages Au atoms tend to migrate 
into the substrate so that a complete Ag layer 
is formed on the top of the Au adsorbate. 
They have also shown that at coverages, B, 
beyond a monolayer (1ML) of Au interdiffu- 
sion no longer happens but Au atoms tend 
to form clusters on the top of the Ag layer 
on the surface. Thus Au grows on Ag (110) 
in a layer+island or Stranski-Krastanov (SK) 
growth mode. The first principles calcula- 
tions by Chan et al. || confirm that it is en- 
ergetically favourable for Au atoms to be cov- 
ered by one Ag layer at submonolayer regime. 
This has been reproduced in Surface Embed- 
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ded Atom Model (SEAM) based Molecular 
dynamics (MD) simulations, which also show 
an onset of SK growth mode for coverages be- 
yond 9 = IML !• 

In this Letter the temperature and cov- 
erage dependence of Au/Pd(110) system de- 
picted from LEED studies of Schmitz et. al 
|T0| is studied theoretically. The structural 
phase transformations show four main fea- 
tures: (i) at sub monolayer coverages only 
(lxl) structures are observed, (ii) for rs 
1.5ML there is a transition from (lxl) struc- 
ture to a mixture of (1 x 2) and (1 x 3) do- 
mains at T rs 310K and back to (1 x 1) at 
T « 53CFC, (iii) For 6 « 2ML a sequence of 
transitions (lxl) ^> (lxn) (1x2) ^» 
(lxl) is seen, where (lxn) denotes a mixture 
of (1 x 2) and (1x3) domains; (iv) increas- 
ing coverage beyond = 2ML leads to ever 
fainter LEED patterns suggesting clustering 
of Au atoms on the top of two Au layers, 
i.e. growth in SK mode. The related LEED 
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studies by Kaukasoina et al. |! l] mainly con- 
firm the previous observations apart from two 
exceptions: they were not able to reproduce 
the (1 x 2) pattern for 6 = 1.5ML and for 
G = 2.0ML they found a continuous tran- 
sition from (1 x 1) pattern via (1 x 2) to a 
mixture of (1 x 2) and (1x3) patterns and 
then back to (1 x 1) structure. In addition, 
for = 2.5ML they found a (1 x 2) recon- 
structed surface with a concentration profile 
(1,0.9,0.36,0.67,0.35), i.e., atomic concen- 
tration of Au is 100% in the first layer, 90% 
in the second, etc. 

A feasible theoretical method to incorpo- 
rate vibrational motion into surface free en- 
ergies is the local Einstein approximation to 
quasiharmonic theory as formulated by Le 
Sar et al. 
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where Ep is the potential energy and \Di\ 
is the determinant of the 3x3 matrix 
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rr^fr- At surfaces, where the three- 
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dimensional symmetry breaks down, it is 
more advantageous to use the local Einstein 
approximation with off-diagonal elements of 
D ia g than a Debye-type approximation [ITS 



which utilizes the trace of D ia/3 and neglects 
shear properties. It has been shown in Ref. 
HH in a study of grain boundaries in Si, that 



although Eq. (1) may fail in the absolute val- 
ues of free energies, the excess free energies 
between two different structures are remark- 
ably good. An optimal structure is obtained 
by allowing the system relax to a local mini- 
mum of Eq. (1). In this study a modification 
of Polak-Ribiere optimization method as ex- 
pressed in Ref. |T5|] , is used with numerical 
approximants of the gradient of free energy 
with respect to atomic positions. 

It has been shown that the Sutton-Chen 
(SC) potential predicts the reconstructions of 
(llO)-surfaces of d-metals correctly except for 
/r |. In addition, SC-potential is easy to 
apply in Eq. (1), in practice, and there ex- 
ist simple but physically satisfactory rules to 
describe bonding between atoms of different 
species. SC-potential |16j is written as fol- 
lows 
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where (m,n) = (7,12) and (8,10), e = 
4.1790 x lO~ 3 eV and 1.2793 x 10~ 2 eV, a = 
3.89A and 4.08A, and c = 108.27 and 34.408 
for Pd and Au, respectively. In heteroge- 
nous y45-systems the following rules m AB = 



AB 



\{n A + n B ), 



,AB 



\{m A + m B ), n 
V a A a B and e AB = V e A e B are found to give 
satisfactory structural, energetic and elastic 
properties fT7|j . 

The calculations were carried out using a 
slab of 14 layers in the [110]-direction, each 
layer consisting of 6 rows and 6 columns 
(36 x 14 = 504 atoms altogether) with peri- 
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odic boundary conditions in the direction of 
both rows and colums. The size of the calcu- 
lational cell is limited by the cut-off radius of 
the SC potential, r c = 2a, and the number of 
rows must be divisible by 2 and 3 for (1 x n) 
reconstructions with n either 2 or 3 and the 
slab must be thick enough to model substrate 
bulk properties deeper within the slab. In 
[110]-direction four layers at the very bottom 
of the slab are fixed in order to maintain bulk 
structure. The atoms are labeled according 
to their species, thus each atom can be de- 
fined to be either of Au or Pd. Rows in top 
layers are removed to create (1x2) and (1x3) 
missing row structures. 

Before optimizing the surface structure, 
the slab is allowed to expand uniformly to 
optimize the bulk free energy of Pd at each 
temperature separately. After this the sur- 
face free energy is minimized by allowing lo- 
cal relaxations using Polak-Ribiere algoritm. 
The surface free energy is calculated by sum- 
ming up the free energies of the atoms in the 
first eight layers subtracted by the bulk free 
energy of Pd per atom and dividing the sum 
by the area of the surface. 

The surface free energies of the follow- 
ing ordered structures were studied: a bare 
(1 x 1) surface of Pd, (1 x 2) reconstructed 
structure with 6 = 0.5ML of Au and an un- 
reconstructed surface with = 1ML of Au. 
For G = 1.5ML an ordered (1 x 2) structure 
with a concentration profile (1,1,0) is stud- 
ied as well as an unreconstructed structure 
with a profile (1, |, 0) with a disordered sec- 
ond layer (and with non-zero configurational 
entropy). In the case of O = 2ML there are 
four simple ordered structures : (lxl) with 
the profile (1,1,0) @ (Fig. 1. (a)), (1 x 2) 
with the profile (1, 1, |) (Fig. 1. (b)) and two 
possibilities for (1x3), either with the profile 
(1,1,1) (Fig. 1. (d)) or (1,1, |, |) (Fig. 1. 
(c)) with Au distributed along the saw-tooth 
shaped pattern forming effectively one atom 
thick layer along the corrugation. The latter 



(1x3) configuration proved more favourable. 
An unreconstructed layer with B = 3ML of 
Au is compared with the corresponding struc- 
ture of O = 2ML to see whether layering 
on the top of the second layer should take 
place. To discuss the concentration profile 
found in Ref. |TT| for a (1 x 2) reconstructed 
film of = 2.5ML two concentration pro- 
files, (1, 1, 1, 0) and (1, 1, 0, 1), are compared. 
In addition, some disordered structures for 
G = 2ML are analyzed to discuss the effect 
of configurational mixing entropy. 

Five main observations to give a frame- 
work for the experimental results can be 
found (See Fig. 2): (i) At submonolayer 
range Pd atoms in the neighbourhood of the 
vacuum have a very high energy. Thus the 
free energies of a bare Pd surface and a (1 x 2) 
surface with O = 0.5ML are quite high in 
free energy compared to (1 x 1) structure of 
a monolayer of Au. Although the bulk free 
energies of Au and Pd are quite close to each 
other, the relatively long range of Pd — Pd 
interactions make it energetically very expen- 
sive to expose Pd atoms to the vacuum. One 
monolayer with (1 x 1) structure, whose sur- 
face free energy is used as the reference in Fig. 
2, is the stablest structure up to room tem- 
perature where (1x2) reconstructed surface 
with = 1.5ML becomes virtually equally 
stable. 

(ii) For = 1.5ML at low temperatures 
a reconstructed (1x2) surface has a lower 
energy than an unreconstructed surface (Fig. 
2.) but, the low energy structure (1 x 1) 
with = 1ML may lead to a possibility 
of (1 x 1) domains at low temperatures. In- 
creasing temperature beyond T « 500K the 
reconstructed and unreconstructed surfaces 
with = 1.5ML have virtually equal surface 
free energy. However, the unreconstructed 
structure should gain in configurational en- 
tropy at high temperatures, since the ran- 
domness of Au atoms in the second layer has 
little cost in free energy. Taking into account 
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this effect, it is to be expected that a LEED 
pattern of (1 x 2) structure for = 1.5ML is 
observed at some narrow temperature range 
(as in Ref. 0). 

(iii) For © = 2 ML there is a sequence 
of transitions as seen in Fig. 2.: (1 x 1) -> 
(1 x n) -> (1 x 2) -> (1 x 3) -> (1 x 1) 
with n being a mixture of 1, 2 and 3, which is 
something very similar to what was observed 
in Refs. PUP] , apart from incorrect temper- 



atures. The phase diagram is complicated 
due to many factors contributing to the ther- 
mal free energy, obviously having very differ- 
ent temperatures of dominance. For exam- 
ple, the higher order the reconstruction is, 
the more corrugated the surface, and increas- 
ing corrugation makes an Au film effectively 
thinner, thus bringing Pd atoms closer to the 
vacuum, which increases surface free energy. 
On the other hand, corrugation decreases the 
proportion of Au atoms at sites with low co- 
ordination, which decreases surface free en- 
ergy. In addition, missing rows alter the lo- 
cal Einstein frequencies, which also has some 
complicated effects on free energy. Some care 
must be taken in interpreting these results, 
since the coverage is thick enough to facili- 
tate interdiffusion. 

(iv) Coverages beyond = 2ML are 
energetically expensive except at high tem- 
peratures, where = 2.5ML is seen to de- 
crease surface free energy (Fig. 2). Although 
not shown in Fig. 2, the unreconstructed 
= 3 ML surface is very high in free en- 
ergy throughout the temperature range and 
the same applies to (1 x 3) reconstructed 
film although the latter is clearly less ex- 
pensive than the unreconstructed surface at 
high temperatures. This result suggests that 
layer-like growth beyond = 2ML does not 
take place and the faint LEED patterns at 



thick coverages |L0| are caused by clustering 
of Au on the top of the first two layers. 

(v) At low temperatures both concen- 
tration profiles, (1, 1, 1, 0) and 



(1x2) reconstructed film with = 2.5ML 
are virtually equally expensive in free energy 
(see Fig. 2). Both the configurations ap- 
proach lower coverages in surface energy as 
the temperature is increased with a devia- 
tion in favour of the profile without an inter- 
mediate Pd layer. However, the difference is 
so small that inclusion of intermixing could 
change the figure. Thus the present analysis 
does not rule out the possibility of interme- 
diate layers with high concentration of Pd at 
high coverages as suggested in Ref. [jTT| . The 
main condition for such an intermediate layer 
to exist is a thick Au coverage that prevents 
Pd atoms from being exposed to the vacuum. 

The main source of inaccuracy here is the 
absence of intermixing, which means ignoring 
configurational mixing entropy. For alloys, it 
is possible to deal with the mixing entropy us- 
ing Bragg- Williams type approximation, as is 
done in Ref. 
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1,1,0,1) for 



In that study the ensemble 
is grand canonical, while here it is canonical 
due to a fixed coverage. But even in this case 
the mixing entropy per layer per unit area for 
randomly mixed layers can be approximated 
as follows: 

S k = -nk B [c k ln(c k ) + (1 - c k )ln{l - c k )} 

(3) 

where c k is the concentration of Au atoms in 
the k th layer and n is the number of atoms 
per unit area. Disorder may increase vibra- 
tional free energy, but the increase may be 
overcome by the contribution of mixing en- 
tropy, E k TS k . 

The effect of mixing entropy is roughly 
estimated for = 2ML by calculating the 
free energy, Eq. (1) and the entropy, Eq. (3), 
for a couple of energetically favourable con- 
centration profiles with disorder. For an un- 
reconstructed surface the profile (1, |, \) has 
a low cost in free energy despite two ran- 
domly mixed layers (The Au atoms of the 
2 nd Au layer are diluted to the 2 nd and the 
3 rd layers). For (1x2) reconstructed sur- 
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faces a corresponding low cost/high entropy by SERC grant No. GR/H58278) of the De- 
profile would be (1,1,5), but with only one partment of Materials in Oxford University 
randomly distributed layer, leading to a half is gratefully acknowledged, 
of the entropy of (1 x 1) surface. For (1x3) 
surface only the atoms in the 4 th layer of 
the profile (1, 1, |, |) are allowed to be re- 
distributed without a high expense, and the 
gain in entropy is slightly smaller than for the 
mixed (1 x 2) structure. In Fig. 3 are shown 
the lowest limits of free energies corrected 
by mixing entropies for mixed layers with for 
randomly distributed Au atoms as described 
above. The complicated transition sequence 
survives the error analysis, because the con- 
figurational entropy affects to the same direc- 
tion for all the structures, although favour- 
ing (1 x 1) surface and disfavouring (1 x 3) 
surface. This might change the sequence of 
transitions more like (lxl)— > (1 x n) - ► 
(1x2)— > (lxl), which actually is what Ref. 
|10| suggests. In any case, the effect of con- 
figurational entropy is significant above room 
temperature. 

Keeping in mind the restrictions of the 
model the present study gives results in qual- 
itative agreement with experimental obser- 
vations. The calculations do not merely re- 
produce the structural phase transitions ob- 
served by LEED but also provide explana- 
tions to some of the experiments. The cal- 
culations take only a couple of minutes per 
point on an HP workstation and thus the 
method can easily be used in planning ex- 
periments and interpreting LEED or STM 
results. Including intermixing and optimiza- 
tion of the concentration profile to the model 
is a logical and necessary continuation to the 
present work. 
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FIGURES 



FIG. 1. Different surface structures for 
9 = 2ML of Au (black circles) on Pd(llO) sur- 
face (open circles): (a) (1 x 1), (b) (1 x 2), and 
(c) and (d) (1 x 3) structures as seen from the 
side of the slab in the direction of rows. The 
structure of (c) has lower free energy than the 
structure of (d). 

FIG. 2. Excess free energies for different cov- 
erages and surface structures, the surface free 
energy of (1 x 1) structure of G = 1ML being 
the reference. For (1x3), the structure of Fig. 1 
(c) has been used. In the legend, 2.5ML (non.) 
refers to the nonmonotonous concentration pro- 
file of G = 2.5ML. 

FIG. 3. The error analysis for excess free en- 
ergies of different structures for 6 = 2ML. The 
symbols without guiding lines show the lower 
limits of the free energy as correted with the con- 
figurational entropy. 
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